Charge density correlations in t-J ladders investigated by the CORE method 
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Using 4-site plaquette or rung basis decomposition, the CORE method is applied to 2-leg and 
4-leg t-J ladders and cylinders. Resulting range-2 effective hamiltonians are studied numerically on 
periodic rings taking full advantage of the translation symmetry as well as the drastic reduction 
of the Hilbert space. We investigate the role of magnetic and fermionic degrees of freedom. Spin 
gaps, pair binding energies and charge correlations are computed and compared to available ED and 
DMRG data for the full Hamiltonian. Strong evidences for shori range diagonal stripe correlations 
are found in periodic 4-leg t-J ladders. 
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Competition between superconducting correlations 
and charge ordering has long been a challenge to nu- 
merical computations ^ of low-dimensional strongly 
correlated electron systems. Spin and hole-doped lad- 
ders 1^ offer an ideal system to investigate the cross-over 
between one to two dimensions. The 2-leg ladder for ex- 
ample is known to exhibit a robust spin gap at and close 
to half-filling as well as hole pair binding |4[ ||]. Domi- 
nant power-law d^2_y2-like pairing and Akp charge den- 
sity wave (CDW) correlations at small doping are charac- 
teristic of a Luther-Emery (LE) liquid regime [|[ Q . How- 
ever, the spin gap magnitude drops sharply as the num- 
ber of legs is increased, e.g. from 0.50 J in the Heisen- 
berg 2-leg ladder to 0.190 J in the Heisenberg 4-leg lad- 
der Both the increase of the magnetic correlation 
length (of the undoped ladder) as well as the drastic re- 
duction of the available ladder length for increasing leg 
number restrict enormously the accuracy of standard nu- 
merical techniques like Exact Diagonalisation (ED) and 
Density Matrix Renormalisation (DMRG) techniques. In 
addition, the DMRG method is limited (in pratice) to 
Open Boundary Condition (OBC) in the leg direction. 

In this Rapid Communication, we use the Contractor 
Renormalisation (CORE) method 0, |l0) to investigate 
hole-doped 2-leg and 4-leg ladders ||l^ . Our aim is to get 
further insights on the issue of pairing and density cor- 
relations from the investigation of large enough systems 
with Periodic Boundary Conditions (PBC) in the ladder 
direction. Such investigations are greatly needed to com- 
plement available DMRG calculations using OBC. Our 
approach is done in two steps; (i) first, using an appro- 
priate partition into small subsystems, we use the CORE 
method to construct an effective hamiltonianx which in- 
tegrates out quantum fluctuations at short length scales; 
(ii) we use ED techniques (supplemented by finite size 
analysis) to compute various physical properties (pair 
binding, spin gaps, etc..) and compare them to those 
of the original model. The method provides new results 
such as e.g. strong evidences for stripes correlations in 
translationally invariant 4-leg ladders. 
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FIG. 1: CORE decomposition in term of plaquette or rung 
sub-systems; (a) 2-leg ladder split in plaquettes; (b) 4-leg 
ladder split in 2 x 2 plaquettes; (c) 4-leg cylinder split in 
4-site rungs. 



We shall consider here a generic n-leg t-J ladder, 

+ *lcg X!(4aC«-Hl,a + h.C.) + irung ^{c\ aCi^a+l + h.C.) , 



where ct^a are projected hole fermionic operators. Open 
(ladders) or closed systems (cylinders) along the rungs 



with isotropic coupling, iieg — trung = 1 and Jiog = 
"Aung — J, will be of interest here. 

In order to implement the CORE algorithm the lad- 
ders are decomposed in small 4-site sub-units as shown 
in Fig. |l] whose M low-energy states are kept to define a 
reduced Hilbert space. The full hamiltonian (||) is then 
diagonahsed on N connected units (with OBC) to re- 
tain its low-energy states. These true eigenstates 
are then projected on the reduced Hilbert space (tenso- 
rial product of the M states of each unit) and Gram- 
Schmidt orthonormalized An effective hamilto- 
nian containing N-body interactions with identical low- 
energy spectrum can then be constructed in terms of the 
reduced basis by a unitary transformation ||ll|] . For sake 
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of simplicity, we shall restrict ourselves to the range-2 
approximation(7V = 2) [T^ . 

At half-filling, retaining in each 4-site unit only the 
lowest singlet and triplet states (4 states) gives excellent 
results Away from half-filling, the simplest trunca- 
tion, referred to as "B" approximation, is to include, in 
addition, the lowest singlet hole pair on the 4-site unit (of 
d-wave symmetry in the case of a plaquette). Formally, 
one can define 4 hard-core bosonic (plaquette or rung) 
operators describing the 4 possible transitions from the 
singlet half-filled GS (vacuum) to one component of the 
triplet state, ta^i, or to the hole pair state, bi (Ref. 13). 



The effective B-hamiltonian TC^ can then be written as 
a sum of a simple bilinear kinetic term Ti'' + 7i* and a 
quartic interaction H"* (Refs. |l^, |l|), 
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where (ij^jOg creates two triplets on plaquettes i and j, 
which are coupled into total spin S. Such an efi'ective 
hamiltonian may serve for analytic and numerical treat- 
ments. Its parameters listed for J — 0.35 and J = 0.5 in 



table p 
model 



are consistent with those found for the Hubbard 



Although, gives already a faithful descrip- 
tion of the physics of the original model, a systematic 
improvement can be done by adding to the above local 
basis 4 extra "fermionic" states corresponding to the de- 
generate {Sz = ±1/2, even and odd chirality or parity) 
single hole GS of the 4-site unit (hereafter referred to as 
"BF" approximation). 

The 2-leg t-J ladder offers an ideal system to test the 
efficiency of the CORE method and the choice of the 
plaquette decomposition. As seen from the behavior of 
the pair-binding energy Ap ~ 2Eo{nh = 1) — i?o("-h = 
2) — Eo{nh = 0) plotted in Fig. ^(a) and from the pla- 
quette charge density-density correlation in the two hole 
GS of the BF- hamiltonian plotted in Fig ^(a), pairs are 
found to be strongly bound and localized almost on a sin- 
gle plaquette. This confirms a posteriori the relevance of 
CORE and of the local basis. Furthermore, the finite size 
scaling of the spin gap for a fixed number of = 2 holes 
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FIG. 2: Finite size scaling analysis vs 1/L for a periodic 
2-leg 2x1/ t-J ladder at J = 0.5 using the decomposition 
of Fig. |l](a) and the effective B- or BF- hamiltonians (as in- 
dicated on plot). DMRG & ED data for the original t-J lad- 
der are also shown for comparison. ED data obtained with 
odd ladder lengths are averaged over boundary conditions (see 
[p^). DMRG data and L ^ oo ED extrapolations are shown 
by arrows, (a) Pair- binding energy Ap. (b) Spin gap of the 
2-hole doped ladder, (c) Spin gap of the 1/8-doped ladder. 




FIG. 3: Correlations as a function of distance (in units of the 
original bond length) up to r = Z//2 in 2 x L t-J ladder at 
J = 0.5. (a) Plaquette charge density correlations for L = 24. 
The BF-hamiltonian (B-hamiltonian) is used for = 2 holes 
(otherwise), (b) Hole pair density- correlation in the lowest 
triplet excited state of 2 x 20 and 2 x 24 2 hole-doped ladders 
using the B-hamiltonian. 
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TABLE I: Parameters of Ti.^ (in units of t) computed for the 
t-J ladder model using range-2 CORE with 2 plaquettes (row 
1 and 2) or 2 4-site rungs (row 3). 




FIG. 4: Finite size scaUng analysis vs 1/L for a periodic 
4-leg 4 X L t-J ladder (open symbols) using the decomposi- 
tion of Fig. §(b) for J = 0.35 (full lines) and J = 0.5 (dashed 
lines). Data are also shown in the case of a periodic cylin- 
der (filled symbols) for J = 0.35 using the decomnosition of 
Fig. blc). Symbols and notations similar to Figs. H(a-c). (a) 
Pair-binding energy Ap. (b) Spin gap of the 2-hole doped 
system, (c) Spin gap of the 1/8-doped system. 



(see Fig. |(b)) or at 1/8 hole density (Fig. |2|(c)) gives 
gaps with 10-20% accuracy in comparison to existing nu- 
merical data Due to the small size of the hole pairs, 
accurate results are obtained even when fermionic exci- 
tations are not included. Note that the effective models 
lead to a smooth finite size behavior, in contrast to the 
original t-J model where "band-filling" effects may lead 
to oscillatory behaviors as seen in Figs. ^ 

We point out the qualitative agreement between our 
results and those of Siller et al. Q who used a more in- 
volved hard-core charged boson model with longer range 
repulsive interactions (giving rise to a Luttinger liquid be- 
havior), but neglected both fermionic and gapped triplet 
excitations [ p^ . Our more systematic and general treat- 
ment using the B-hamiltonian gives a similar qualitative 
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FIG. 5: Hole-pair density- density correlation on a 4 x 12 lad- 
der at J/t = 0.35. PBC are used in the leg direction and 
correlations are measured from the reference plaquette on the 
lower left corner. From top to bottom, = 4, 6, 8. The 
surfaces of the dots are proportional to the values of the cor- 
relations. 



picture as can be seen from the charge correlations shown 
in Fig|^(a); we observe the characteristic 4ki?-CDW spa- 
tial oscillations of the LE phase showing the same number 
of maxima as the number of hole pairs. Let us empha- 
size that this is also in agreement with DMRG calcu- 
lations [T^ . Our approach performed on finite homo- 
geneous systems is then complementary to the DMRG 
technique using OBC. 

Although the agreement with the hard-core charged 
boson model is qualitatively good, we believe that includ- 
ing magnetic triplet excitations in the local basis is never- 
theless important to describe interplay between magnetic 
and pairing correlations. For example, it is known that 
the lowest triplet excitation in a 2 hole doped (or very 
weakly doped) t-J ladder consists of a hole pair-magnon 
boundstate Indeed, the extrapolated value of the 

spin gap in the presence of 2 holes (see Fig. ^(b)) is lower 
than that of the undoped ladder (0.5 J) and than the hole 
pair binding energy (shown in Fig. ^(a)). Moreover, as 
seen in Fig.||(b), the correlation between the hole pair 
density and the plaquette S^-component clearly shows 
an enhancement at short distance fl^ . 

We finish the investigation of the 2-leg ladder by using 
the effective Hamiltonian to calculate the Luttinger liq- 
uid parameter Kp which governs the long distance power- 
law behavior of the charge correlations related to the 
unique massless charge mode. Some values of Kp ob- 
tained from the Drudc weight D and the compressibility 
K [|l5j as Kp — TT^y Dk/2 are listed in table In ad- 
dition, we also list here the charge velocity Up obtained 
from the relation Up = nD/ Kp and which agrees within 
a few percents to the values obtained directly from the 
linear dispersion of the charge mode. Note also that 
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these values compare very well to existing ED [Q and 
DMRG [ITtI data. 
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TABLE II: Parameters Kp and function of dop- 

ing computed on 2 x 24 (") and 2 x 28 (") ladder with B- 
hamiltonian and J /t — 0.5. 

We now turn to the investigation of the 4-leg t-J lad- 
der (with OBC along rungs) or cylinder (with PBC along 
rungs) , for which the best choices of unit decomposition 
are depicted in Fig. |l|(b) and Fig. |^(c) respectively. Re- 
sults for pair binding energies and spin gaps are shown 
in Fig. ^(a-c). Results for ladders and cylinders are sim- 
ilar although the hole pair binding is much stronger in 
cylinders where hole pairs are preferably formed on cross- 
sectional plaquettes (periodic rungs) rather than on "sur- 
face" plaquettes. Generically we found that the pair 
binding energy is larger than the spin gap of the un- 
doped (Heisenberg) system (0.190J for the 4x L ladder). 
Therefore, the lowest triplet state in the 2-hole doped (or 
very lightly doped) 4-leg ladder is similar to a Heisenberg 
ladder magnon, which may be (or may not be) loosely 
bound to a hole pair depending whether its excitation 
energy is lower or equals the magnon energy of the un- 
doped system. Since the data shown on Fig. |^(b) are not 
fully conclusive we have computed in addition the hole 
pair density-S'2 correlation and found, as for the 2-leg 
ladder case, an enhancement of the spin density on the 
neighboring sites of the hole pair suggesting, indeed, the 
existence of a hole pair-magnon boundstate. 

Upon increasing doping, as seen from the hole-pair 
density-density correlation shown in Fig.|^, we observe a 
clear tendency of the hole pairs to align along the diago- 
nal (1, ±1) directions with a periodicity corresponding to 
one pair every two plaquettes, a behavior also reported in 
DMRG calculations [p^ |l9| and reminiscent of the picture 
of diagonal stripes. Note that real space charge correla- 
tions are fully consistent with the power law decay found 



in the effective charge boson model |19 . 

To conclude, the CORE method is a powerful method 
to extract effective hamiltonians for strongly correlated 
models. It allows numerical simulations on significantly 
larger systems than those available for the original model. 
We show that including charge and spin bosonic excita- 
tions gives reliable results as long as the hole pair binding 
energy is not too small. Results for the effective model 
of the 2-leg t-J ladder are in excellent agreement with 
known analytic and numerical data. Within the effective 
models hole pair-triplet bound states form for both 2-leg 
and 4-leg ladders, a key feature to be compared to 5*0(5) 
phenomenological theories |20|. In addition, the method 



enables unbiased (since calculated on translationnaly in- 
variant clusters) analysis of hole pair density correlations. 
While 4fcF-CDW correlations are found in 2-leg ladders, 
our computations provide clear evidences in favor of short 
range diagonal stripes in 4-leg ladders. 
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